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Abstract 

We study the non-boost-invariant evolution of a quark-gluon plasma subject to large early- 
time momentum-space anisotropies. Rather than using the canonical hydrodynamical ex- 
pansion of the distribution function around an isotropic equilibrium state, we expand around 
a state which is anisotropic in momentum space and parameterize this state in terms of three 
proper-time and spatial-rapidity dependent parameters. Deviations from the Bjorken scal- 
ing solutions are naturally taken into account by the time evolution of the spatial-rapidity 
dependence of the anisotropic ansatz. As a result, we obtain three coupled partial differ- 
ential equations for the momentum-space anisotropy, the typical momentum of the degrees 
of freedom, and the longitudinal flow. Within this framework (0+l)-dimensional Bjorken 
expansion is obtained as an asymptotic limit. Finally, we make quantitative comparisons 
of the temporal and spatial-rapidity evolution of the dynamical parameters and resulting 
pressure anisotropy in both the strong and weak coupling limits. 

Keywords: Non-Boost-Invariant Dynamics, Anisotropic Plasma, Non-equilibrium 
Evolution, Viscous Hydrodynamics. 



1. Introduction 

In recent years relativistic viscous hydrodynamical models have been applied with great 
success to the description of the non-central anisotropic flow measured at the Relativistic 
Heavy Ion Collider llr|8|. Despite this success, it is well known viscous hydrodynamics is 
not an accurate description during early times after the initial impact. The hot and dense 
matter created after the collision is rather small in transverse extent and expands very 
rapidly in the longitudinal direction. In addition, transverse expansion of the matter is not 
expected to become important until times on the order of 3 or 4 fm/c after the collision. As a 
consequence, large momentum-space anisotropies are developed at early-times. When these 
anisotropies are large the shear viscous corrections can be of the same order or larger than 
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the isotropic pressure, casting doubt on the validity of viscous hydrodynamics to faithfully 
model the dynamics. For example, it is known that within the framework of 2nd-order 
viscous hydrodynamics large momentum-space anisotropies can lead to negative longitudinal 
pressure at early times j^j. 

One is therefore motivated to develop an alternative approach that will allow us to better 
understand the transition between early-time anisotropic expansion and late-time hydrody- 
namical behaviour in ultrarelativistic heavy ion collisions. The work presented here fills this 
gap using a technique which allows one to smoothly connect both regimes using a unified 
theoretical framework. We present a method which allows us to derive hydro dynamical- like 
evolution equations in the presence of large momentum-space anisotropies. The method con- 
sists of changing the usual hydrodynamic expansion of the one-particle distribution function 
in the local rest frame 

/(x,p,r) = / eq (|p|,T(r)) + 5/ 1 + 5/ 2 + ... , (1) 

which is an expansion around an isotropic equilibrium state / eq (|p|,T), to one in which the 
expansion point itself can contain momentum-space anisotropies 

/(X, P, r) = /aniso(P,Phard, 6) + 8f{ + & f'2 H ■ ( 2 ) 

In Eq. (j2J) £ is a parameter that measures the amount of momentum-space anisotropy and 
Phard is a non-equilibrium momentum scale which can be identified with the temperature of 
the system only in the limit of isotropic equilibrium. We introduce an ansatz for the leading 
order term of Eq. fl2]), /aniso(p,Phard> £)> which approximates the system via spheroidal equal 
occupation number surfaces as opposed to the spherical ones associated with the isotropic 
equilibrium distribution function [10j. We do not calculate the remaining terms of the 
corrections 5f' n , their evaluation is beyond the scope of this paper. Nevertheless, it may be 
possible to perform such a calculation in a consistent way for relativistic system with the 



use of irreducible anisotropic tensors [ll| or spheroidal harmonics |12J for which our ansatz 
is the leading order term. Moreover, one expects that the corrections 8f' n will have smaller 
magnitude than the isotropic corrections Sf n because of the momentum-space anisotropy is 
built into the leading order term in the reorganized expansion. 

In this paper we extend our previous work on highly anisotropic dynamics [13) to include 
both temporal and spatial-rapidity dependence of the system thereby removing the assump- 
tion of boost invariance of the system. By making use of an ansatz for /aniso(p,Phard> 0> 
we derive leading order evolution equations for a non-boost-invariant plasma which expands 
only in the longitudinal direction. Based on this ansatz we obtain three coupled partial 
differential equations for the momentum-space anisotropy, the typical momentum scale of 
the particles and the longitudinal flow rapidity. We show that these equations reduce to 
the boost-invariant description of the plasma in the appropriate limit. Finally, we present 
numerical solutions of the non-boost-invariant evolution equations and make a quantitative 
description of the pressure anisotropy for weak and strong coupling. 
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2. Kinetic theory approach for a non-boost-invariant and highly anisotropic 
QGP 



In this section we derive the equations of motion necessary to describe the transition be- 
tween early time dynamics and late-time hydrodynamical behaviour of a non-boost-invariant 
system. For simplicity, we assume that there is no transverse expansion and that all expan- 
sion is along the longitudinal direction]]] The derivation is based on taking moments of the 
Boltzmann equation. Before presenting the evolution equations, we first review the moment 
method of the Boltzmann equation applied to general relativistic theories 15Ml7||. then we 
apply this formalism to the case of a non-boost-invariant system which is anisotropic in 
momentum space. 

2.1. Moments of the Boltzmann equation 

For a dilute system, the Boltzmann relativistic transport equation for the on-shell phase 



space density /(x, p, t) is 15 



p"0„/(x,p,t) = -C[/] 



(3) 



where C[f\ is the collisional kernel which is a functional of the distribution function. In 
order to have a tractable approach to dissipative dynamics from kinetic theory, one usually 
takes moments of the Boltzmann equation ([3]) to determine the equations of motion of the 
dissipative currents [151-117 . 

The equation of motion for the n-th moment of the Boltzmann equation ia^| 



o Tfli(12... Hn+i 



where 



Jp 

f p^ p ^ ... p ^c[f] . 
Jp 



(4) 

(5a) 
(5b) 



Eq. (|3]) defines an infinite set of coupled equations for these moments. By truncating the 
expansion, one obtains a finite set of equations fl5j |. In the Israel-Stewart (IS) theory, 
the evolution equation for the dissipative currents and the transport coefficients can be 
extracted from the second moment of the Boltzmann equation within the 14 Grad's ansatz 
for the distribution function 16, 17 . 



1 We point out that the method presented here can be generalized to include azimuthal and non-azimuthal 
transverse expansion. We postpone these studies to a future publication [l4| • 
2 Our notation for the momentum-space integral is J = (27r) -3 J d 3 p/p°. 
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By using Eq. (jlj) for n = 0, 1, 2; one finds 



/ I?dpf(x,p)=d lt N'>= [ C[f] 



(6a) 



p -'p 

p^d,f(x,p) = d ll T tlv = [ p"C[f] = , (6b) 



/ p^p%f(x, P ) EE ^J^ A = / P VC[/] EE P^ A . (6C) 
•/ p «/ p 

The zeroth moment of the Boltzmann equation (IS"a|) . tells us the evolution of the particle 
current N^. When f p C[f] = 0, we have particle number conservation, however, for gen- 
eral gauge theories like QCD where inelastic processes are important, away from chemical 
equilibrium we have J f C[f] ^ 0, i.e. there is a source term for particle production and 



annihilation [18142 ll|. The first moment of the Boltzmann equation (16bl) . gives us the con- 
servation of the energy and momentum since f p^Clf] = 0. Eq. (|6cl) represents the balance 
of the dissipative fluxes. I^ vX is a symmetric tensor of the dissipative fluxes and P uX is a 



traceless tensor that is related to dissipative processes due to collisions [l5|, \Y7\, |22|. The 
right hand side of Eq. (I6"c|) does not trivially vanish, so in general f p u p x C[f] = P uX ^ 0. 

To obtain the complete set of equations for the dynamical variables such as the particle 
density n, energy density S, etc., in hydrodynamical treatments one canonically expands 
the distribution function around equilibrium by using some functional form for the fluctua- 
tions (see Eq. ([!])). In our approach, we do not expand the distribution function around an 
isotropic expansion point. Instead, we find evolution equations for the kinematical param- 
eters contained in the leading term / an i SO in Eq. by taking moments of the Boltzmann 
equation. 

2.2. Energy-momentum tensor for a highly anisotropic plasma 

The general energy-momentum tensor for an anisotropic plasma is given by BH 



T nu = ( E + Vt) u h u * _ Vt g n» _ (<p T _ r L )v^v u , (7) 

where £, Vt, and Vl are the energy density, transverse pressure, and longitudinal pressure. 

is a unitary four-vector which defines the longitudinal flow velocity of the system. In 
the local rest frame (LRF), u^ 1 = (1,0,0,0). In addition to the fluid velocity, we introduce 
a vector that is a normalized space-like vector (t> M t> M = — 1) which is orthogonal to 
{u^ = 0). In the LRF, = (0, 0, 0, 1)1 

The energy momentum tensor takes a diagonal form T^ v = diag(£, Vt, Vt, Vl) in the 
LRF. The evolution of its components are given by energy-momentum conservation 

dfjT^ = . (8) 



3 Usually in hydrodynamics, one uses operator A Mi/ = g^ v — u^u u which is orthogonal to the four-velocity 
vf*. From its definition, it is straightforward to see that is an eigenvector of the operator A Miy with 
eigenvalue equal to one, i.e. t^ iV v v — . 



If we project the last expression with and v M we obtain 24 



VS = - (S + V T )6 + (V T - Vl^Vv" , (9a) 
pp L = _ (7? T _ v L ) e+(S + Vt) , (9b) 

where V = w M <9 M , V = v^d^, 6 = and 6 = d^. 

2. 3. Boltzmann equation for the anisotropic distribution function 

In order to make analytical progress, we will assume that the leading term of Eq. (j2j) is 
given by the ansatz proposed by Romatschke and Strickland [10]. Within this approach the 



leading order anisotropic distribution function in the LRF can be obtained from an arbitrary 
isotropic distribution function (/i SO ) by squeezing or stretching / iso along one direction in 
momentum space 

/(x, p, r) = f 1S o(p^r(x)/p 2 hard (x)) , (10) 

where Phard is related to the average momentum of the partons and ^ v (x) is a symmetric 
tensor which indicates the direction and strength of the anisotropy. This particular choice 



is a generalization of the original longitudinal RS ansatz 10 . 

In this work we restrict ourselves to a system which expands along the longitudinal 
direction in a non-boost-invariant way so that ^ v takes the form 

r»=diag(l,0,0,£(x)), (11) 

where — 1 < £ < oo is a parameter that reflects the strength and type of anisotropy. Using 
this, the first argument of Eq. (flOl becomes 

r{x)p»p v = (po) 2 +& 2 



z > 



= Pt+(1 + 0p:, (12) 

where we have considered on-shell massless particles and p\ = p 2 . + p 2 . Choosing ^ v as 
given by dTTj) allows us to recover the original form of the RS ansatz as used by the authors 
of Ref. [lji.e. 

/(X, P , r) = / RS (p, £,Phard) = / iso ([PT + (1 + 0^]/Phard) ■ ( 13 ) 

To write the non-boost-invariant Boltzmann equation for the RS ansatz ffTB"]) . it is convenient 
to switch to the coordinates specified by r = y/t 2 — z 2 and tanhg = z/t. In this coordinate 
system, the RS ansatz ([13]) is 

/its(p,£,2W,tf) = /iso(p^[l + (1 +£)smh 2 (y-tf)]/pLd) , (14) 



4 The anisotropy parameter £ is related to the average longitudinal and transverse momentum of the 
plasma partons via the relation 2|| 26 1 £, = ^(Pt)/(Pl) ~ 1- The system is locally isotropic when £ = 0. 
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where i? is the hyperbolic angle associated with the comoving or LRF and y is the lab 
momentum-space rapidity. Note that £, phard and i? are functions of both r and <;. After 
changing variables, the Boltzmann equation is 

p T cosh(y - + PT Smh ( y - ^ £j / R s(p,e,Phaxd^) = -C[/ RS ] . (15) 

So far we have not prescribed a specific collisional kernel C[f}. Since the calculation of the 
collisional kernel is difficult in general, we model it here by considering the relaxation time 
approximation 

CN=p/r [/Rs(p,£,Phard,tf)-/ eq (|p|,T)] , (16) 

where / eq (|p|,T) is the local equilibrium distribution and T is the relaxation rate. The tem- 
perature T is determined dynamically by requiring energy conservation £ eq (r) = £non-eq( T ) Ui 



27] . From the results of Appendix A one finds for the RS ansatz ffT3l) that T = -R 1//4 (£)Phard [l 



281 ] . As we demonstrated in a previous paper [L3|], T can be related to the shear viscosity 
and shear relaxation time of the system by matching our treatment to 2nd-order viscous 
hydrodynamics in the small anisotropy limit. We will return to this point in the end of the 
next section. 

2.4- Evolution equations for the kinematical parameters of the RS ansatz 

In this section we apply the moment method to the anisotropic RS ansatz (fT3j) for a 
non-boost-invariant system. The boost invariant case was already studied in a previous 



publication [13J . We will assume in what follows that the system is conformal and therefore 



2. 4-1- Oth moment of the non-boost invariant Boltzmann equation 

As we mentioned in Sect. I2.1[ the Oth moment of the Boltzmann equation is equivalent 
to the equation for the evolution of the particle current (Eq. (J6"a|) ). After plugging the 
RS ansatz (fT51) into the Boltzmann equation, the left hand side (LHS) of Eq. (|15p is 

LHS = dJ\M\^{Ed T ^ + ^d^)-^4^{ 1 +^( Ed ^ + -d^) 

L Kard K T J p 2 h&vd V T J 

— [E <9 T p hard + — <9 f Phard ) , (17) 

where oj = \p\ + (1 + C)Pz) /Phard aric ^ we h ave rewritten the result in terms of the energy 
and longitudinal momentum in the LRF defined in the fluid cell via E = px cosh(?/ — $) 
and p z = pTsmh(y — respectively. Taking the zeroth moment (integrating over the 
momentum-space measure (27r) -3 J d 3 p/E) of the last expression together with the right 
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hand side of Eq. (fl5|) . we obtain 



(p 2 z d w f iso (u)) 



T 



Phavd 



E 1 T Hard 
^ CVPhard 



2(up z d UJ f iso (u)) 

Phard 



ard 



-r 



Phard 

(/iso(p,f,Phard)) - </ eq (p,T(r))) 



where (a) = (27r) 3 J d 3 p. These integrals can be calculated analytically using the general 
form of the RS ansatz (TT3"j) 



(Pa^/isoM) 



1 Phard 



2(l + 3/2 



™iso(P 



hard J 







{Ep z duf- lso {u))) = 0, 



(w0 u /tao(w)) 

(/Rs(P,^,Phard)) 



3 n iso (p h 



ard J 



2 (1 + 1/2 

o, 

^-iso (Phard) 



(i + 1/2 ' 

</ cq ( P ,T(r) = ^(OPhard)) = ^ 3/4 (0^iso(Phard) 



(19a) 

(19b) 
(19c) 

(19d) 
(19e) 
(19f) 
(19g) 



where in the last line we have used the Landau matching condition T(r) = 7£ 1//4 (£)phard- 
After replacing these averages in Eq. (1181) . we obtain our first partial differential equation 



6 



Phard 



-d T Ph. 



ard 



2T 



(20) 



2.4-2. 1st moment of the non-boost invariant Boltzmann equation 

In Sect. 12.11 we pointed out the equivalence between the first moment of the Boltzmann 
equation and energy-momentum conservation (Eq. f!6bl) ). Therefore, we can use Eqs. fl9]) 
derived from the energy- momentum tensor given by ([7]). For the longitudinal dimension, 
it is useful to use the following parametrizations for the vectors and introduced in 

sect. uHSi 



= (cosh0(r,?),O,O,sinh0(T,?)) 
= (sinhtf(T,0, 0,0, coshi?(r,?)) 



(21a) 
(21b) 
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where i? is the hyperbolic angle associated with the velocity of the LRF as measured in the 
lab frame. With this particular choice, one finds 

u^Vv" = 9= (sinh(tf - q) d T + cosh ^ ~ d ^ # ; (22) 

Vfl Vu^ = -6 = - (cosh(?9 - q) 8 T + Smh ^ - ^ d,y . (23) 

Using these identities, in this coordinate system we can rewrite Eqs. ([9]) as^| 

(ft+ tMhW ;- t) - t) »)g(r, t )= 

- (5(r, + P L (r, ? )) (tanh(tf(r, q) - q) d T + -£J0(t, q) , (24a) 
(tanh(tf(T, 5 ) - C )S T + -f)v L (T,l) = 

- (£(r, + PUr, 0) (» + ta " h(tf ^ ~ B.) < (r, <) ■ (24b) 

Using the RS ansatz (1131) it is possible to obtain analytic expressions for the energy density 
and longitudinal pressure (see Appendix A Eqs. OA.laP and 0A.lc| . respectively). Using 
these results, plugging them into Eqs. (|24p and performing some algebra, we finally obtain 
the remaining two evolution equations 

fc'(0 q t + 4 ^rPhard + tanll(?9 - q) ( K' (£) g + 4 l\p harc A 
^(0 T Phard T f Phard / 

= - (l + ) ( tanh (^ -S)dr + ^y, (25a) 

_ (3 m +l)( ^^ as) ,, (25b) 

where and 7Zl(0 are defined in Eqs. f lA.laj) and ( lA.lcj) . respectively, while Tl'(q~) and 
TZ' L (0 are 



n' L (0 = 



If 1 — £ arctan\/£\ 

^ (1 + 2- £3/2 J 

l + £) arctan\/£ £ — 1 



2 £3/2 2£(l + £)" 

Eqs. (J22D together with (|2"01) give us three coupled partial differential equations for the 
parameters £(r, q), i?(r, <;) and Phard( T , of the RS ansatz (TT3"j) . Their evolution will describe 



5 A similar structure of these equations was first found by Satarov et al. 29|, 3pj ] . 
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the transition between early-time dynamics and late-time hydrodynamical behaviour for a 
non-boost-invariant plasma in a unified framework. 

The (0+l)-dimensional boost invariant case can be recovered within our approach if we 
require that i?(r, q) = q, £(t, <j) = £(t) and Pha,rd( T , <?) = Phard( r )- By taking this limit in 
Eqs. (1201) and f[2l)l) and rearranging some terms, we obtain two differential equations for £ (r) 
and PhardO) 



1 



6 



no 



Phaid 

4 



-d T p h 



ard 



d r i~ H <9 r p h 

Phard 



ard 



2r 
i 

T 



i-^ 3 / 4 (07TTe' 



I 1-1 

£(1 + 0^(0 £ 



(26a) 
(26b) 



The last expressions correspond precisely to the boost invariant case studied before by us 
(see Eqs. (20) and (22) in Ref. Il3|). For the boost invariant case, one can also show that 
free streaming and ideal hydrodynamical behaviour are reproduced when T = and r — > oo, 
respectively 13]. In addition, the (0+l)-dimensional 2nd-order viscous hydrodynamics of 
Israel and Stewart can be recovered from our ansatz in the small £ limit if we identify 



T = — 



5 T) 



(27a) 
(27b) 



We will implement this matching in our numerical simulations of the non-boost-invariant 
evolution equations ( 12"U|) and (l25p. In this manner the evolution of the system is specified 
by the shear viscosity with no need to introduce additional parameters, e.g. thermalization 
time scales, etc. 



3. Results and Discussion 

In this section we present the results of numerically integrating the coupled nonlinear 
differential equations using the relation between the relaxation rate T and the ratio of shear 
viscosity to equilibrium entropy fj = T]/S which results from Eqs. ( 12"T|) . namely 



r = 2T(r) = 27^/ 4 (Op hard 
which should be used for T in Eqs. fj20|) and f l25|) . 



(28) 



3.1. Initial Conditions 

In order to minimize the number of figures presented here we will only consider LHC-like 
initial conditions. Additionally, for purposes of this paper we will present two types of initial 



6 We refer the reader to Sect. 2.4 of Ref. [l3[ for further details. 
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conditions: (I) an initially isotropic distribution which has a Gaussian profile in the number 
density and (II) an initially anisotropic distribution which has the same Gaussian density 
profile (I). For the number density profile in spatial rapidity (<j) we use a Gaussian 

which successfully describes experimentally observed pion rapidity spectra from AGS to 
RHIC energies |3ll435j| and extrapolate this result to LHC energies. The parametrization we 



use is 



with 



= n exp 



2(j2 



3 (1 - c*) 



In (^/s NN /2m p ) 



(29) 



(30) 



where c s is the sound velocity, m p = 0.938 GeV is the proton mass, for LHC ^s NN = 5.5 TeV 
is the nucleon-nucleon center-of-mass energy, and Uq is the number density at central rapidity. 
In this paper we will use an ideal equation of state for which c s = 1/ v3 in natural units. 

In case (I) we can use Eq. f l29|) to straightforwardly calculate the initial dependence of 
Phard on rapidity by using Eq. flA.2j) with r = r ) = to obtain 



Phard(?,7- = 0) 




(31) 
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Figure 2: Evolution of the dynamical variables in the case of isotropic Gaussian initial conditions and strong 
coupling viscosity corresponding to fj = l/(47r). All panels show evolution at central rapidity, except panel 
showing time evolution of the hyperbolic angle 9 — ^ which shows the evolution at c = 5. 



where po = Phard(? = 0,r = 0) is the initial "temperature" at central rapidity. For LHC 
conditions we use po = 845 MeV at an initial time of tq = 0.1 fm/c. In case (II) we choose 
£(5, r = t ) = 100 which means that we should increase the central value of Phard by a factor 
of (1 + Q 1 ' 6 — 2.16 in order to start with the same initial particle density distribution [see 
Eq. (1A.2j) ]. Doing this gives in case (II) p Q ~ 1.82 GeV. 

In both cases (I) and (II) we must also fix the initial condition for the hyperbolic angle 
$ which specifies the four- velocity of the local rest frame. Here we choose $(<;, r = tq) = <; 
following Refs. 0, 0, M ■ This is not the most general initial condition possible, of course; 



however, absent some guiding principle we choose here to follow Refs. [24J, [29|, [30[ and 
approximate the initial state as having a four-velocity profile corresponding to an initially 
boost-invariant configuration. 

3.2. Case I: Isotropic Initial Conditions 

In this subsection we present results for the spatial-rapidity and proper-time dependence 
of our dynamic parameters in the case that the parton distribution function is assumed to 
be locally isotropic at r = tq = 0.1 fm/c. We present the cases of transport coefficients 
corresponding to a strongly-coupled plasma and a weakly-coupled plasma separately 
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Figure 3: Evolution of the rapidity dependence of the dynamical variables in the case of isotropic Gaussian 
initial conditions and weak coupling viscosity corresponding to fj = 10/(47r). 



3.2.1. Strongly Coupled 

In Fig. [T] we plot the spatial-rapidity dependence of the anisotropy parameter (upper 
left), the ratio of the longitudinal to transverse pressure (upper right), the hard momentum 
scale (lower left), and the hyperbolic angle of the local rest frame four- velocity as measured 
in the lab frame (lower right). We have fixed the ratio of the shear viscosity to entropy 
density to be fj = rj/S = l/(47r). Each line in the plot corresponds to a fixed proper- 
time t G {0.1, 2.1, 4.1, 6.1} fm/c. From these figures we see that large momentum-space 
anisotropies are developed at large rapidity. This can be traced back to the fact that the 
initial temperature is lower at large rapidity, causing the rate of relaxation back to isotropy 
to be slower @. At r = 6.1 fm/c we see that in the case of strong-coupling shear viscosity 
at q = 5, V L jV T ^ 0.75. 

We note that the deviation of °d ^ q which results for r > r is indicative of the breaking 
of longitudinal boost invariance. The fact that d > <j results from there being pressure 
gradients in the longitudinal direction which increase the longitudinal velocity beyond that 
which would be obtained if the flow were boost invariant. The fact that the •& — q curves turn 
over at large spatial rapidity is an indication that the material at the edges is propagating 
at a slower longitudinal velocity than the material in the central region. We also note that 
as shown in the lower left plot of the hard momentum scale both the width and amplitude 
of the momentum distribution are changing with proper-time. The is width increasing with 
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Figure 4: Evolution of the dynamical variables in the case of isotropic Gaussian initial conditions and weak 
coupling viscosity corresponding to fj — 10/(47r). All panels show evolution at central rapidity, except panel 
showing time evolution of the hyperbolic angle 9 — ^ which shows the evolution at c = 5. 



proper-time due to longitudinal pressure gradients. 

In Fig. [2] we plot the proper-time dependence of our dynamical variables and the pressure 
anisotropy for fixed spatial rapidity. We show the anisotropy parameter at central rapidity 
(upper left), the ratio of the longitudinal to transverse pressure at central rapidity (upper 
right), the hard momentum scale at central rapidity (lower left), and the hyperbolic angle 
of the local rest frame four- velocity as measured in the lab frame at 5 = 5 (lower right). We 



do not show $ at central rapidity since $(^ = 0) = by symmetry 29|, |30j. As one can see 



from the plots presented in Fig. [2] at central rapidities large momentum-space anisotropies 
can develop at early times even in the case of a strongly-coupled plasma. From the upper 
right panel of Fig. [2] we see that at r ~ 0.25 fm/c one has Vl/Vt ^ 0.5. Such large pressure 
anisotropies could have a significant effect on observables which are sensitive to the early 



time dynamics of the quark-gluon plasma such as photon [361-142 
and heavy-quarkonium production 451454 



dilepton |25|, |2J, |43|, |44 



3.2.2. Weakly Coupled 

In Fig. [3] we plot the spatial-rapidity dependence of the anisotropy parameter (upper 
left), the ratio of the longitudinal to transverse pressure (upper right), the hard momentum 
scale (lower left), and the hyperbolic angle of the local rest frame four- velocity as measured 
in the lab frame (lower right). We have fixed the ratio of the shear viscosity to entropy 
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Figure 5: Evolution of the rapidity dependence of the dynamical variables in the case of anisotropic Gaussian 
initial conditions and strong coupling viscosity corresponding to fj — l/(4.7r). Initial value of £ is £o = 100, 
independent of rapidity. 



density to be fj = rj/S — 10/ (Ait) . Each line in the plot corresponds to a fixed proper- 
time t G {0.1, 2.1, 4.1, 6.1} fm/c. From these figures we see that large momentum-space 
anisotropics are developed at all rapidities. At r = 6.1 fm/c we see that in the case of weak- 
coupling shear viscosity Vl/Vt ^ 0.4 at all rapidities. In terms of the hyperbolic angle 
we see weaker deviations from boost-invariant longitudinal flow than seen in the case of 
strong-coupling transport coefficients (Fig. [1]). Such strong momentum-space anisotropies 
indicate that a naive viscous hydrodynamical treatment would not be trustworthy at any 
spatial rapidity shown. This is in line with the results presented in Ref. H where we 
discussed the generation of negative longitudinal pressure due to large momentum-space 
anisotropies when using 2nd-order viscous hydrodynamical evolution equations. Using the 
partial-differential equations presented here we are able to evolve the system even in the 
case of large momentum-space anisotropies and the longitudinal pressure is guaranteed to 
be positive at all times. 

In Fig. H] we plot the proper-time dependence of our dynamical variables and the pressure 
anisotropy for fixed spatial rapidity. We show the anisotropy parameter at central rapidity 
(upper left), the ratio of the longitudinal to transverse pressure at central rapidity (upper 
right), the hard momentum scale at central rapidity (lower left), and the hyperbolic angle 
of the local rest frame four- velocity as measured in the lab frame at q = 5 (lower right). 
As one can see from the plots presented in Fig. H] at central rapidities large momentum- 
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Figure 6: Evolution of the dynamical variables in the case of anisotropic Gaussian initial conditions and 
strong coupling viscosity corresponding to fj — l/(47r). All panels show evolution at central rapidity, except 
panel showing time evolution of the hyperbolic angle 9 — <; which shows the evolution at <j = 5. Initial value 
of £ is £o = 100, independent of rapidity. 



space anisotropies persist at all times. From the upper right panel of Fig. H] we see that 
at r ~ 0.25 fm/c one has Vl/Vt ^ 0.25. This is much larger than the strongly-coupled 
case presented in the previous section. This suggests that anisotropic photon, dilepton, 
and heavy-quarkonium production rates integrated in spacetime using our dynamical model 
could be used to constrain quark-gluon transport coefficients. This could be done by com- 
paring theoretical predictions resulting from the assumption of strong-coupling (Fig. [1]) and 
weak-coupling (Fig. |3]) transport coefficients and integrating over the resulting spacetime 
evolution. Finally we point out that one sees from the lower left panel of Fig. H] that the 
evolution of the hard momentum scale is far from that expected due to weakly-viscous 
expansion. Initially we see a period of approximate free-streaming expansion which later 
transitions to viscous hydrodynamical expansion. This would have the effect of increasing 
the lifetime of the quark-gluon plasma in the weakly-coupled case and could result in stronger 
transverse flow than naively predicted by a 2nd-order viscous hydrodynamical treatment. 

3.3. Case II: Anisotropic Initial Conditions 

In this subsection we present results for the spatial-rapidity and proper-time dependence 
of our dynamical parameters in the case that the parton distribution function is assumed 
to be locally anisotropic at r = tq = 0.1 fm/c. As detailed in the beginning of the results 
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Figure 7: Evolution of the rapidity dependence of the dynamical variables in the case of anisotropic Gaussian 
initial conditions and weak coupling viscosity corresponding to fj = 10/(47r). Initial value of £ is £o = 100, 
independent of rapidity. 



section we take £(r ) = 100 which corresponds to a highly oblate distribution with V L Vt- 
Such highly oblate distributions are predicte d by color-glass-condensate models of the early- 



time dynamics of a heavy-ion collision [19|, |55|. As in the previous section that assumed 



isotropic initial conditions, we present the cases of transport coefficients corresponding to a 
strongly-coupled plasma and a weakly-coupled plasma separately. 

3.3.1. Strongly Coupled 

In Fig. we plot the spatial-rapidity dependence of the anisotropy parameter (upper 
left), the ratio of the longitudinal to transverse pressure (upper right), the hard momentum 
scale (lower left), and the hyperbolic angle of the local rest frame four- velocity as measured 
in the lab frame (lower right). We have fixed the ratio of the shear viscosity to entropy 
density to be fj = rj/S = 1/(4-71"). Each line in the plot corresponds to a fixed proper- 
time t G {0.1, 2.1, 4.1, 6.1} fm/c. From these figures we see that large momentum-space 
anisotropies persist at large rapidity At r = 6.1 fm/c we see that in the case of strong- 
coupling shear viscosity at <j = 5, Vl/Vt — 0.7. This is similar in magnitude to the result 
obtained when one starts with an isotropic initial condition (see Fig. [1]). We also note a 
characteristic triangular shape to the spatial-rapidity dependence of the pressure anisotropy 
(Fig. 0top right). 

In Fig. Owe plot the proper-time dependence of our dynamical variables and the pressure 
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Figure 8: Evolution of the dynamical variables in the case of anisotropic Gaussian initial conditions and 
weak coupling viscosity corresponding to fj = 10/(47r). All panels show evolution at central rapidity, except 
panel showing time evolution of the hyperbolic angle 9 — <; which shows the evolution at ? = 5. Initial value 
of £ is £o = 100, independent of rapidity. 

anisotropy for fixed spatial-rapidity. We show the anisotropy parameter at central rapidity 
(upper left), the ratio of the longitudinal to transverse pressure at central rapidity (upper 
right), the hard momentum scale at central rapidity (lower left), and the hyperbolic angle 
of the local rest frame four- velocity as measured in the lab frame at <j = (lower right). 
As one can see from the plots presented in Fig. at central rapidities large momentum- 
space anisotropies persist in the case of a strongly-coupled plasma. From the upper right 
panel of Fig. |6]we see that at r ~ 0.25 fm/c one has Vl/Vt ^ 0.5. This is similar to the 
case of isotropic initial conditions suggesting a kind of "attractor" for the momentum-space 
anisotropy evolution. This attractor can be identified as the Navier-Stokes solution if one 
analyzes the late time asymptotic solutions of the evolution equations presented here and/or 
2nd-order viscous hydrodynamics. 

3.3.2. Weakly Coupled 

In Fig. [7] we plot the spatial-rapidity dependence of the anisotropy parameter (upper 
left), the ratio of the longitudinal to transverse pressure (upper right), the hard momentum 
scale (lower left), and the hyperbolic angle of the local rest frame four- velocity as measured 
in the lab frame (lower right). We have fixed the ratio of the shear viscosity to entropy 
density to be fj = rj/S = 10/(47r). Each line in the plot corresponds to a fixed proper- 
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time r G {0.1, 2.1, 4.1, 6.1} fm/c. From these figures we see that large momentum-space 
anisotropics are developed at all rapidities. At r = 6.1 fm/c we see that in the case of 
weak-coupling shear viscosity at q = 5, Vl/Vt — 0.2. 

In Fig. [8] we plot the proper-time dependence of our dynamical variables and the pressure 
anisotropy for fixed spatial rapidity. We show the anisotropy parameter at central rapidity 
(upper left), the ratio of the longitudinal to transverse pressure at central rapidity (upper 
right), the hard momentum scale at central rapidity (lower left), and the hyperbolic angle 
of the local rest frame four- velocity as measured in the lab frame at q = (lower right). As 
one can see from the plots presented in Fig. [8] at central rapidities large momentum-space 
anisotropies persist for a long time in the case of an initially anisotropic weakly-coupled 
quark-gluon plasma. From the upper right panel of Fig. [S] we see that at r ~ 0.25 fm/c one 
has Vl/Vt ^ 0.01. Such extreme momentum-space anisotropies would preclude the use of 
a viscous hydrodynamical treatment; however, our reorganized expansion offers some hope 
to treat systems with such high anisotropies. 

4. Conclusions and Outlook 

In this paper we have derived three coupled partial differential equations given in Eqs. (12"U|) 
and fl25l) whose solution gives the time evolution of the plasma anisotropy parameter £, the 
typical hard momentum scale of the partons phard; and the longitudinal flow velocity vari- 
able We relaxed the assumption of boost-invariance by allowing the longitudinal flow 
velocity variable to differ from the spatial rapidity q and we also allowed the hard momen- 
tum scale and plasma anisotropy to depend on the spatial rapidity. The partial differential 
equations were obtained by taking moments of the Boltzmann equation using a relaxation 
time approximation collisional kernel and a spheroidally anisotropic ansatz for the under- 
lying partonic distribution function. By requiring that our equations reduced to 2nd order 
viscous hydrodynamics in the limit of small anisotropy we were able to obtain an analytic 
connection between the relaxation rate T appearing in the collisional kernel and the plasma 
shear viscosity r\ and shear relaxation time r^. Using Eqs. fl20l) and fl25l) and T as a functions 
of r) we were able to numerically solve the coupled partial differential equations in both the 
strong and weak coupling limits. 

Our numerical results indicate that in both the strongly and weakly coupled cases large 
momentum-space anisotropies can be generated. Therefore, a treatment such as the one de- 
rived here in which momentum-space anisotropies are built into the leading order ansatz can 
better describe the time evolution of the system. Relatedly we demonstrated that within our 
reorganized approach at all spatial rapidity the longitudinal pressure remains positive during 
the entire evolution. This is to be contrasted to 2nd-order viscous hydrodynamics which 
can give negative pressures indicating the breakdown of the expansion about an isotropic 
equilibrium state. We have also discussed the fact that our equations can reproduce both 
extreme limits of the dynamics: ideal hydrodynamical expansion when the shear viscosity 
goes to zero and free streaming when the shear viscosity goes to infinity. 

The numerical results obtained here differ from those obtained in Ref. [24j in that we see 
a much slower relaxation to isotropy at large spatial rapidity. This results from the fact that 
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the authors of Ref. 24j assumed that there is a single rapidity-independent thermalization 
time scale which they arbitrarily chose to be 0.25 fm/c. In this work we instead implemented 
a method for matching to 2nd-order viscous hydrodynamics that was introduced by us in 



Ref. |13|. The conclusion of that paper was that the relaxation rate T should be proportional 
to hard momentum scale (see Eq. (I28p herein). If in a certain region the average momentum 
scale is lower, as is the case in the forward regions described by large spatial rapidity this 
implies a slower relaxation to isotropic equilibrium. In our case this naturally arises due to 
the matching to 2nd-order viscous hydrodynamics in the limit of small anisotropy and as a 
result it not necessary to introduce external parameters such as thermalization times as a 
function of rapidity. 

Beside the interesting theoretical goal of improving the description of systems which have 
large momentum-space anisotropies, we point out that the dynamical equations derived 
here can be used to assess the impact of rapidity-dependent moment-space anisotropies 
on observables which are sensitive to the early-time dynamics of the quark-gluon plasma 



such as photon [36H42| . dilepton [25l. l26l. l43l |44I| and heavy- quarkonium production [45H54 



Additionally, one could consider using the evolution of £ and Phard presented here in order 
to fold in a realistic anisotropy evolution into simulations of non-abelian plasma instabilities 



1Q]|56|-|63|. 



Looking forward it is of critical importance to include the transverse expansion and ellip- 
tical flow of the system. This can be done minimally by allowing for one more momentum- 
space anisotropy parameter that describes the momentum-space anisotropy in the transverse 
plane and allowing all dynamical parameters to depend on proper-time, spatial rapidity, and 
transverse position. Work along these lines is currently underway 
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Appendix A. Analytic expressions for the components of the energy momen- 
tum tensor, number density, and entropy density 



The energy-momentum tensor T^ v = (2tt) 3 J rf 3 p/p° p^jf f(x , p,t) for the RS ansatz 
(fT3"|) is diagonal in the comoving frame and its components are 13j, |28| 

If 1 arctanv^A 
£{Phavd,0 = T ^2\T+~C + ~= J ^iso(Phard) , (A.la) 

= 7£(£) ^iso(Phard) , 

Pt&w, = 1 (t xx + ryy) = A + j M0 ) n so (p har d) , (A.ib) 

= ^ T (on so (Pha rd ) , 

V L ( Phald , = -I? = | p + n so (phard) ; (A lc) 

= ^ L (On S °(Phard) , 

where "P^°(phard) and PjfG^hard) are the isotropic transverse and longitudinal pressures and 
£iso(Phard) is the isotropic energy density. 

One can also easily obtain the number density as a function of phard and £ 

^,P h a rd ) = ^^OC^f=, (A.2) 

V 1 + <5 V 1 + £ 

where n; so is the isotropic number density. Similarly, one can obtain for the entropy density 

of/- \ ^iso(Phard) / » \ 

S (£, Phard) = ■ ( A ' 3 ) 
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